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Abstract. The dynamics of dislocation assemblies in deforming crystals indicate the 
emergence of collective phenomena, intermittent fluctuations and strain avalanches. 
In polycrystalline materials, the understanding of plastic deformation mechanisms 
depends on grasping the role of grain boundaries on dislocation motion. Here the 
interaction of dislocations and elastic, low angle grain boundaries is studied in the 
framework of a discrete dislocation representation. We allow grain boundaries to 
deform under the effect of dislocation stress fields and compare the effect of such 
a perturbation to the case of rigid grain boudaries. We are able to determine, 
both analytically and numerically, corrections to dislocation stress fields acting on 
neighboring grains, as mediated by grain boundary deformation. Finally, we discuss 
conclusions and consequences for the avalanche statistics, as observed in polycrystalline 
samples. 



PACS numbers: 62.20.Mk,05.40.-a, 81.40Np 



Dislocation interactions mediated by grain boundaries 



2 



1. Introduction 

Recent advances in experimental techniques have allowed the redefinition of the once 
supposedly well-understood problem of crystal plasticity. In spite of the traditional 
picture of a smooth fluid-like motion, acoustic emission (AE) measurements on 
deforming ice crystals [U [2] as well as more recent compression tests on metallic 
specimens [3j 0] have during the last decade made it clear that plastic deformation 
of crystalline materials proceeds as a heterogeneous and intermittent adjustment to 
applied external drives. Such intermittence in plastic deformation had been for long 
observed [51 E] and dismissed as considered unworthy of serious consideration. The 
recent discoveries challenge this picture since that heterogeneity appears to take over 
through the emergence of a burst-like scale-free behavior. Plastic deformation is found 
to develop through strain avalanches, characterized by power law distributed sizes and 
energy emissions. Slip motion is found to follow fractal patterns [7] and deformed 
specimens exhibit self-affine surfaces [8]. 

The observation of scale invariance suggested that fluctuations should not be 
imputed to the discrete nature of the deforming lattices. Instead, the collective behavior 
of crystal dislocation was held responsible for the emergence of scale-free patterns, in a 
close-to-critical behavior fashion. Such an observation triggered a number of numerical 
[21 El EE]] and, whenever feasible, analytical [HI [12] studies so that the validity and the 
limits of applicability of this picture could be assessed. 

At present, there is evidence that the notion of scale invariance in crystal plasticity 
as a consequence of a close-to-critical state is consistent with observation in single and 
multi-slip geometries, over a variety of materials, hardening coefficients and loading 
conditions. However, the general understanding of these processes is still chiefly 
phenomenological and several aspects of such phenomenology are still obscure. In 
particular, it is still not entirely clear to which extent such considerations should be 
extended to the case of polycrystalline geometries. 

Grain structure has always been supposed to play a central role in crystal yielding 
[T3| [H] , and has attracted considerable attention also more recently [T5], [161 CEO HH1 HH] • 
The so called Hall-Petch law, which relates the yield stress of a polycrystal to its average 
grain size, has been recently reformulated in the light of collective dislocation dynamics 
models [20J and its breakdown for nanometer sized grains explained as a result of the 
failure of those systems to initiate or sustain such a process. 

These observations pose the general question on how strain avalanches spread 
through grain boundaries (GB). The most simplistic approach consists in assuming that 
the elastic interactions of two dislocations on two sides of a GB are entirely screened by 
the GB. Such a crude assumption holds only if the different slip geometries of neighboring 
grains allow that. In general, however, nothing prevents dislocation stresses to act on 
neighboring grains, at least in a dampened form due to geometrical incompatibilities. 

A reasonably balanced approach consists in postulating that plastic flow 
(i.e. dislocation dynamics) within a certain grain builds up back stresses due to 
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incompatibility of orientations and slip configurations with the neighboring grains. As 
soon as the polycrystal reaches the yield point, the accumulated back stresses are 
redistributed within the neighboring grains. Plastic yield is thus accompanied by a 
relaxation process by which dislocations effectively interact through grain boundaries 
by transferring incompatibility stresses. 

The advantage of such an approach is that its results indirectly account for 
experimental observation. By looking at cumulative probability distributions of emission 
amplitudes in AE experiments [9], for example, creeping ice crystals exhibit a power- 
law decay with a rather universal exponent (3. Polycrystalline samples still retain the 
power-law behavior, however with a lower (and possibly non-universal) exponent and, 
as expected, a large size cut-off due to finite grain sizes. A natural explanation of the 
cut-off is that to some extent strain avalanches are confined within grains. However, 
explorations of avalanche models [9] seem to suggest that the decrease of the exponent 
P is possible only if excess stresses within one grain can eventually trigger dislocation 
avalanches in the neighboring ones. 

Dislocation interactions through grain boundaries are thus essential to explain 
plastic flow in polycrystals. The idea of incompatibility stresses acting on neighboring 
grains certainly captures part of the physics of the problem but it leaves several questions 
unanswered, by virtue of its generality. For instance, the exact mechanism by which 
stresses are relaxed is not clear, as well as not much can be said about the typical range 
of stress rearrangements on a quantitative basis. 

In this paper, we perform an analytical and numerical study of dislocation 
interactions through grain boundaries, in the light of the above considerations. In 
order to make the problem treatable, we consider a relatively simple GB model, that 
of a low-angle grain boundary, or dislocation wall, in a two-dimensional geometry. This 
serves as a starting point to capture the relevant physics of the problem and expose 
it in a conveniently simple way. In passing, we also recall that dislocation walls are 
frequently encountered in deforming crystals as well as in a variety of similar systems 
such as colloidal crystals [21] and flux-line lattices in Type II superconductors [22\. 

In the following we discuss how stresses generated by dislocations deform low- 
angle grain boundaries and how this affects in return the GB stress field perceived by 
dislocations in their proximity. Theoretical predictions are compared with numerical 
results. We find that this perturbative effect accounts for a screening effect for stress 
deforming the grain boundaries. The typical strength of this screening proves inversely 
proportional to the grain size and the perturbation itself is found to be short-ranged 
and exponentially suppressed beyond a distance that again scales with the grain size. 
The analytical and numerical computations are in the next section, while Section [3] is 
devoted to discussion and conclusions. 
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2. GB-dislocation interaction 

Let us consider a small-angle grain boundary of length L in two dimensions. The GB is 
schematized as a discrete linear assembly of N + 1 edge dislocations distributed along 
the y direction, with Burgers vectors b parallel to the x axis, of sign s = ±1. The 
dislocation spacing is assumed constant and set equal to D. The GB is pinned at the 
boundaries (y = and y = L) and can be deformed under an applied shear stress, see 

Fig. m 

By representing the GB displacement through the vector rj, where rj n is the 
displacement of the n— th dislocation along the x axis, the elastic energy cost of a 
given deformation 77 reads 

m - f e (i) 

where = /x6 2 /(47rD 2 (l — z/)). In our notations, is the shear modulus and v the 

Poisson ratio. Sums in Equation (JTJ) run over all couples of dislocations. As a result 

of long range dislocation interactions, Equation ([T]) has the form of a non-local elastic 

functional. Its elastic kernel scales as \k\ in Fourier space, where k is the wave vector 

associated with the GB deformation profile. 

If an external stress field r(r) is applied, forces acting on the GB can be calculated 

in terms of Peach-Koehler interactions. The GB begins to deform and settles in a 

configuration where the applied forces equate the restoring elastic tensions deriving 

form the term [1] above. Given the particular geometry of our problem, the GB profile 

is given by the displacement field i] which satisfies 

oE , . . . . 

j- = rfr(r), (2) 

where the restoring force on the right-hand side is calculated as the (variational in the 
continuum case) derivative of the elastic energy functional in (Tj[|). Carrying on the above 
calculation, Equation (j2J) can be rewritten as a simple linear non-homogeneous problem, 
in the form 



by defining 



^V nm r] m = sbr(r), (3) 



, — K-, — L_ n^m 

V nm = { («-™r ' ^ (4) 

K J2k^n ( n _ fe )2, n = m. 



2.1. A simple case - uniform stress 

Let us first review the simple case in which the applied stress is uniform, i.e. r(r) = r 
|23j. In order to solve Equation ([3]), which can be rewritten in operatorial terms as 

VH = sbr H = V- 1 sbr 1 (5) 
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one needs to diagonalize the matrix of interactions V nm . This corresponds to solving 
the eigenvalue problem 

E Vnmlln = ( 6 ) 



m 

J 



where we have called A/ each of the N eigenvalues and y n each component of the 
eigenvector corresponding to A/. It can be shown that for large enough N, the eigenvalues 
read 

7T 2 l [ib 2 id 1 

Xl = -N2n(l-v)D* =2K N (X L (7) 
and the eigenvectors 

7i = /f sin ^,J = (0), 1,2,3, ...N. (8) 

Such an eigenvector representation is possible since the functions 7^ match boundary 
conditions and constitute a complete set, such that any continuous function with 
piecewise continuous and differentiable derivative can be expanded in terms of them 
|24j . Orthogonality of eigenvectors is ensured by the following condition 

2 . Txl in . nl 2 n 

— > sin sin = 0i,i o - " 

N ~q N N hh K ' 

Now V = T T Ar and V~ x = rA _1 r T , where we have called T the matrix of eigenvectors 
and A the (diagonal) matrix of eigenvalues and we have exploited the orthogonality of 
the transformation T, such that = T T . Easily, we derive H = &7T T A -1 T in the form 

Vn = Sbrj2 KWnlL (10) 
m,l 

and replacing the sum over m with an integral (this is done also in the following - see 
the Appendix for a discussion) 

%\-vDL v- 1 - (-1)' • ?rZn 

r in = s— r T l^ — w — sm ^T' 11 

ti z fi b j 2l z i\ 
where the [1 — (—l) l ]/2 factor simply limits the sum to odd values of the I index. 



2.2. Stress generated by a dislocation 

Solving this problem in the case of a more complicated form of the external stress can 
be a difficult task, instead. In particular, we are interested in the shear stress generated 
by an external edge dislocation a, placed at r a = (x a ,y a ). The dislocation a has a 
Burgers vector of modulus b and sign s a = ±1. Even assuming that the position of such 
a dislocation is fixed, the magnitude of the GB-dislocation interaction will necessarily 
depend on the deformation of the GB itself. Equation becomes 

v nmVm = C n [r a ) x (r) n - x a ), (12) 
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with C n (r a ) = ss a c x [x 2 a - (Dn - y a ) 2 }/[x 2 a + (Dn - y a ) 2 } 2 and c = fib 2 /[2n(l - v)\. 
Here we have assumed x n <C x a , as the non-local stiffness of the GB would not allow 
huge deformations to take place. The equation above can be rewritten in the form of 
Equation as 

X)[Km - C n (r a )]r] m = -x a C n (r a ) (13) 

m 

or, equivalently, 

(V-C)H = -x a C, (14) 

where C is an operator defined accordingly. The obvious solution would be 

H = -x a (V-C)- 1 C, (15) 

however, unlike the case of a uniform stress, the inversion of the (V — C)~ l matrix is 
a tough task. We cannot proceed through a straightforward diagonalization as before, 
since the matrices V and C are never simultaneously diagonal. A feasible way to carry on 
the calculation requires an approximation. Matrix C includes corrections coming from 
the deformation of the GB, which is supposed to be a relatively small effect compared to 
the main interaction carried by V. These considerations can be easily reformulated in a 
more quantitative fashion. As a consequence, using simple relations of matrix algebra, 
we can write (V - C)- 1 = (1 - V~ 1 C)- 1 V- 1 « (1 + V^QV' 1 , which is a first order 
expansion in C . The displacement is then 

H = -a; (r T A- 1 r + T T A^TCT 7, A'^C (16) 

or, in the other notation, 



C m {r a ). (17) 



E KWrrfn + EE Ar 1 A; 1 7 i 7 ^ fe (r a ) 7 h^ 

mil is 

In Equation (IT7|) . the second term within square brackets is considerably smaller 
than the first and can be neglected. The deformation of the GB becomes 

Vn ~ -X a J2K 1 lnT,^mCm(r a ). (18) 
I m 

So far we have made basically two assumptions: i) we deal with a large GB; ii) the 
external dislocation is not in the immediate vicinity of the GB. These choices have 
allowed us to take continuum limits in sums and to disregard higher order terms. At 
this point, however, the sum over m in Equation (|T8l) cannot be calculated easily, unless 
other approximations are made. We restrict ourselves to the limit in which the external 
dislocation sees a quasi-infinite GB, i.e. L ^> x a . The forces exerted by the external 
dislocation on the GB dislocations can be either attractive or repulsive depending on the 
distance and the angle between the external dislocation and each GB dislocation. As a 
result the concavity of the GB profile may change along the y direction. The sum over 
m in Equation ( |T8l) can be replaced by an integral from — oo to +00 (see Appendix). 
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The integral is solved calculating residues around the poles of the C m {r a ) function. The 
displacement thus determined is 

1 ^ _2i| r I Ttlya ixln 
r) n = -2ss a x a — 2_^e i' 0| sin — any. (19) 

This result is compared with numerical results from discrete dislocation dynamics 
simulations of the GB-dislocation system. These are done with a similar setup as e.g. 
in Ref. [25J, with a straight GB as an initial configurations. By letting the GB relax 
in the stress field of an external dislocation, a good agreement is found (see Fig. [2]). 
Notice that in order for the simulation setup to correspond to Eq. (fl9l) (which was 
derived by computing the integral from — oo to +oo), we have considered a system 
where in addition to the deformable part of the GB of length L, there is a large number 
(here 500) of pinned dislocations on both sides of the deformable part of the GB. The 
agreement with Eq. ([191) is the better the larger the number of these pinned dislocations. 

2. 3. Correction to the GB field 

We are finally able to determine how a dislocation induced deformation affects the stress 
field generated by the GB. We consider the case in which L 3> x a , since it matches are 
simulation setups. The shear stress at position (x, y) due to the deformed GB is a 
functional of the displacement field r\ n and can be expressed, by calling a n the shear 
stress of the n— th GB dislocation, as 

^(x,y) = J2 a n(x - r] n ,y - y n ), (20) 

n 

where r\ n is given by Equation [19]. Since we are interested in plastic flow activation 
well inside each grain, we shall assume r] n <C x, that is we keep far enough from the 
stress singularities of GB dislocations. This way we can express stress fields through 
their Taylor expansion around r/ n ~ as o~ n {x — r] n ,y — y n ) w o~ n (x, y — y n ) + cr' n (x, y — 
y n )Vn + ]-/2o~n( x >y ~~ Un)Vn + ^(Vn)- Equation ( 1201) will be rewritten as 

<r(x,y) = aW(x,y) + <r^\x,y) + a^(x,y) + ... (21) 

where the i-th term on the right hand side is obtained summing the i-th order as in 
Equation ( |20l) . In passing we note that, given the nature of the derivation operation, 
in Equation (|2ip even terms will be odd functions of x, while odd terms will be even 
functions of x. Then it is no coincidence that o~(°>(x,y) is the well known stress field 
generated by a straight GB (see References [26J, [27] ) . 

As for the first order o~^(x, y) calculations are not different from the ones carried 
on in the previous sections. We eventually obtain, in the continuum limit of the sum, 

o- {1) (x,y) = Tcs a x a Y^^J2 e ^sin^px (22) 

vriM iily 
e L sin — - — . 



4-* 
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One can observe that this term is, as expected, an even function of x. It accounts for 
corrections observed along the x-direction, where the GB deformation produces an even 
decrease in the absolute value of the shear stress, see Fig. [3] for corresponding results 
from dislocation dynamics simulations. This is balanced by an increase along the 45- 
degree directions along which shear stress is redistributed between adjacent grains. Such 
corrections are introduced by the o~( 2 \x, y) term, which can be calculated integrating the 
second-order term of the expansion of o n [x — i] n ,y — y n ). The procedure does not differ 
from what we have seen so far, with the exception of some additional complications 
arising from the fact that this time the sum over / is squared. Skipping the details of 
the calculation, we obtain 

«*>(*, V) = 2«V.J!L£ £ £ *> ^ - ^ X (23) 

■ , \ \n an2 zdi+MkJ (n \l + Al . . \ ir(l + X)y 

x sign(x) {I + A) e £ — — - — \x\ — 1 cos 

V z ±j J Lj 

„ 7r|t-A||!t| / 7r \l — Al , . \ 7r(Z — X)y 

where, as usual, I, A = 1, 2,3...N. Indeed, this term proves to be an odd function of 
the variable x. The same method allows the calculation of further orders of the stress 
expansion (12T1) which become relevant as test dislocations are close to the GB, where 
the anisotropy of stress fields becomes relevant. 

We should point out that the terms of the sums over I (and A) are exponentially 
suppressed. As a consequence, sums can be stopped at the first terms. In particular, 
if the asymmetry or the exact shape of the GB profile are not a concern, the very first 
terms provide a fair description of the deformed GB. Furthermore, one should notice that 
the correction itself decays exponentially with the distance of the external dislocation 
a. This suggests the emergence of a screening distance, beyond which the effects of 
GB deformation do not affect dislocation-GB interaction anymore. This distance is 
expected to be proportional to L/(irl), i.e. higher order corrections in I are weaker at 
large distance, but may come into play on a shorter range. 



3. Discussion 



Both analytical calculations and numerical studies of the stress fields of a GB-dislocation 
system show that the GB deformation tends to screen the dominating stress that is 
causing the deformation. At the same time, due to the fact that a stress field of an 
edge dislocation has both positive and negative regions, some parts of such a stress field 
are also enhanced in the same process. One should note, however, that the strongest 
effect the GB deformation has is screening close to the deformed part of the GB. This 
happens on both sides of the GB (see Fig. [3]). 

The implications for collective dynamics in polycrystals with flexible grain 
boundaries would then be the following: Low angle GB's tend to accommodate the 
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excess stresses within a grain by deforming in such a way that such excess stresses 
are partly screened. This screening effect is exponentially suppressed beyond a length 
that is of the order of the typical grain size ~ L, and has a magnitude scaling as 
~ 1/L. Consequently, in addition to acting as boundaries to the dislocations themselves, 
deformable GB's tend to hinder also stresses from being transferred from one grain to 
another. This effect gets stronger with decreasing grain size. One should note, however, 
that also various dislocation arrangements present within grains tend to screen the long 
range stress fields typical to these systems. 

The behavior of this simple model system is hence controlled mainly by a single 
parameter, namely the typical grain size ~ L. In fact, our results suggest that 
deformable low angle GB's tend to confine avalanches of plastic activity to some extent 
within single grains. This could be anticipated to be the case also for general (high 
angle) grain boundaries: even if they cannot be thought as deformable arrays of mobile 
dislocations, the large difference in the respective orientations of neighboring grains 
means that only a fraction of the stress is effectively transferred over such a GB. 

Recent experiments on polycrystalline ice reveal an intriguing crossover in the AE 
amplitude distributions: Small and large avalanches appear to be characterized by 
different power law exponents [9]. The observation that small avalanches are described 
by the same exponent as avalanches in single crystal samples leads to the natural 
conclusion that such small avalanches do not feel the presence of the GB's. On the 
other hand, larger avalanches will necessarily interact with the GB's. In the experiment 
presented in Ref. [9], the exponent describing such large avalanches was found to be 
smaller as compared to the one describing the statistics of smaller avalanches. 

Thus, the presence of GB's will change the statistics of avalanches in some way. 
One possibility, along similar ideas as presented in Ref. [9] , is to consider the spreading 
process on a coarse-grained scale of grains. An initial avalanche within a grain, if 
reaching the region close enough to a grain boundary, may spread to neighboring 
grains thus triggering avalanches there. Triggering is governed by an effective spreading 
probability p that depends (among other things) on the details of the dislocation- 
grain boundary interaction discussed above. Hence the ensuing spreading process is 
characterized by an average p and geometrical factors (correlations in GB orientations, 
grain sizes and their correlations etc.). This leads to a number iV of grains that 
participate in the "total" avalanche. The statistics of N will follow at the simplest a 
mean-field branching process for which one in general has three possible scenarios: i) a 
sub-critical process with (N) = const, ii) an exactly critical process with P(N) ~ iV~ 3 / 2 , 
Hi) a supercritical process which never stops. Clearly, the last case can be excluded by 
common sense. As for the the first two options, it appears much more likely that 
the grain-to-grain spreading is actually sub-critical. This would be in agreement with 
the observation that in experiments of Ref. [9] a grain-size dependent cut-off for the 
acoustic emission amplitudes was observed. On the other hand, this would indicate 
that the observed avalanche distributions would not change from what one sees on a 
single grain level. 
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In conclusion, we have considered the GB mediated dislocation interactions for 
a "toy" system in two dimensions. This has allowed us to make analytical progress, 
to compare with numerics, and to understand the nature of screening. Regardless of 
the obvious limitations of the simple model studied here, working out these problems 
leads to an increased understanding of the physics of avalanches in the plasticity of 
polycrystalline materials. Clearly, there is much room for short-term investigations in 
three separate directions at least. One is the accumulation of further experimental data 
on these phenomena. The second is the influence of the geometry (dimension, GB angle, 
grain size, correlations) on these processes, which can be studied partly experimentally 
and partly by numerical, atomistic, simulations. Finally, there seems to be room for 
the development of coarse-grained models which would match their predictions with 
experimental data. 
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Appendix 

In several cases in the derivations, sums have been replaced by integrals to infinity. 
While the validity of the infinity assumption is quite intuitive if L is taken significantly 
larger than x a , the transition to continuum requires a few words of discussion. Provided 
that the sum is extended to infinity, we calculate it using the Poisson integral relation 



Calculating the integral first, we always obtain an exponential ~ e mkx / D term, which 

rapidly goes to zero if k is increased, as we always assumed x 3> D. As a result, the 

k = term of the sum over k alone provides a good estimate of the initial sum. But 

this is equivalent, as it can be easily noticed, to replacing the sum with an integral in 

the first place. 
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Figure 1. A schematic diagram of the grain boundary dislocation system. The GB 
is modelled as a discrete linear assembly of TV + 1 edge dislocations distributed along 
the y direction, with Burgers vectors b parallel to the x axis. Dislocation spacing is 
D. The endpoint dislocations of the GB, positioned at y — and y — L, are fixed to 
x = 0. The rest of the GB dislocations are allowed to move in the x direction as a 
response to the stress field of the external dislocation positioned at x — x a , y = y a . 
The ensuing GB displacement in the x direction is represented by the vector ff, with 
components r\ n corresponding to the displacement of the nth dislocation. 
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Figure 2. Three examples of deformed grain boundary profiles (symbols) in the case 
L = 100 ^> x a . The deformation is due to the stress held of an external dislocation 
positioned at x = x a , y = y a . The solid lines correspond to the predictions of Eq. (1191) . 
500 pinned extra dislocations (with the same spacing as in the deformable part) have 
been inserted on both sides of the GB. Notice that there are no fitting parameters. 




Figure 3. Stress fields of a system consisting of a grain boundary (the dark 
vertical structure) and an external dislocation positioned to the right side of the 
grain boundary (x a — 10, y a — 50, L = 100). Left: Stress field in a situation 
where the grain boundary is not allowed to deform. Middle: Stress field with a 
deformable grain boundary. Right: The difference of the absolute values of the two 
stress fields, |(j de / ormed | _ \ a unde formed^ Negative values indicate screening, positive 
values enhancement of the magnitude of the stress due to grain boundary deformation. 



